Quintessence arising from exponential potentials 
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We demonstrate how the properties of the attractor solutions of exponential potentials can lead to 
models of quintessence with the currently observationally favored equation of state. Moreover, we 
show that these properties hold for a wide range of initial conditions and for natural values of model 
parameters. 
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I. INTRODUCTION 

Measurements of the redshift-luminosity distance re- 
lation using high redshift type la supernovae combined 
with cosmic microwave background (CMB) and galaxy 
clusters data appear to suggest that the present Uni- 
verse is flat and undergoing a period of A driven inflation, 
with the energy density split into two main contributions, 
^matter ~ 1/3 and Q\ ~ 2/3 ffl-S. Such a startling 
finding has naturally led theorists to propose explana- 
tions for such a phenomenon. One such possibility that 
has attracted a great deal of attention is the suggestion 
that a minimally coupled homogeneous scalar field Q (the 
"quintessence" field), slowly rolling down its potential, 
could provide the dominant contribution to the energy 
density today thanks to the special form of the potential 
P,ra. Non- minimally coupled models have also been in- 
vestigated H [n]] . The advantage of considering a more 
general component that evolves in time so as to dominate 
the energy density today, as opposed to simply insert- 
ing the familiar cosmological constant is that the latter 
would require a term pa ~ 10~ 47 GeV to be present at 
all epochs, a rather small value when compared to typical 
particle physics scales. On the other hand, quintessence 
models possess attractor solutions which allow for a wide 
range of initial conditions, all of which can correspond 
to the same energy density today simply by tuning one 
overall multiplicative parameter in the potential. 

There is a long history to the study of scalar field cos- 
mology especially related to time varying cosmological 
constants. Some of the most influential early work is 
to be found in Rcfs. [jjj] [L4J. One particular case which 
at first sight appears promising is the one involving ex- 
ponential potentials of the form V ex exp(A/s<2), where 
k 2 = 8ttG |l2|-|il|. These have two possible late-time 
attractors in the presence of a barotropic fluid: a scal- 
ing regime where the scalar field mimics the dynamics 
of the background fluid present, with a constant ratio 
between both energy densities, or a solution dominated 
by the scalar field. The former regime cannot explain 
the observed values for the cosmological parameters dis- 
cussed above; basically it does not allow for an acceler- 
ating expansion in the presence of a matter background 



fluid. However, the latter regime does not provide a fea- 
sible scenario either, as there is a tight constraint on the 
allowed magnitude of £Iq at nucleosynthesis [ |l7|Jl8| . It 
turns out that it must satisfy Og(lMeV) < 0.13. On the 
other hand, we must allow time for formation of structure 
before the Universe starts accelerating. For this scenario 
to be possible we would have to fine tune the initial value 
of pq, but this is precisely the kind of thing we want to 
avoid. 

A number of authors have proposed potentials which 
will lead to A dominance today. The initial sugges- 
tion was an inverse power law potential ( "tracker type" ) 
V ex Q~ a [p| Jl^Jl9( |, which can be found in models of 
supersymmetric QCD |2yj2l[. Here the ratio of energy 
densities is no longer a constant but pq scales slower 
than pb (the background energy density) and will even- 
tually dominate. This epoch can be set conveniently 
to be today by tuning the value of only one parame- 
ter in the potential. However, although appealing, these 
models suffer in that their predicted equation of state 
wq = pq/ pq is marginally compatible with the favored 
values emerging from observations using SNIa and CMB 
measurements, considering a flat universe [p2|~p4[. For 
example, at the 2<r confidence level in the Qm — wq plane, 
the data prefer wq < —0.6 with a favored cosmological 
constant wq = — 1 (see e.g. |24[]), whereas the values per- 
mitted by these tracker potentials (without fine-tuning) 
have wq > —0.7 [B5[. For an interpretation of the data 
which allows for wq < — 1 see Ref. |2tj| . 

Since this initial proposal, a number of authors have 
made suggestions as to the form the quintessence poten- 
tial could take ]27J-|33|] . In particular, Brax and Mar- 
tin |28] constructed a simple positive scalar potential 
motivated from supergravity models, V oc exp((5 2 )/Q a , 
and showed that even with the condition a > 11, the 
equation of state could be pushed to wq w —0.82, for 
£Iq = 0.7. A different approach was followed by the au- 
thors of f3(j,[33|] . They investigated a class of scalar field 
potentials where the quintessence field scales through an 
exponential regime until it gets trapped in a minimum 
with a non-zero vacuum energy, leading to a period of de 
Sitter inflation with wq — ► — 1. 

In this Brief Report we investigate a simple class of 
potentials which lead to striking results. Despite previ- 



ous claims, exponential potentials by themselves are a 
promising fundamental tool to build quintessence poten- 
tials. In particular, we show that potentials consisting of 
sums of exponential terms can easily deliver acceptable 
models of quintessence in close agreement with observa- 
tions for natural values of parameters. 



II. MODEL 

We first recall some of the results presented in 
||14] , [l7|Jl8| . Consider the dynamics of a scalar field Q, 
with an exponential potential V ex exp(AK<5). The field is 
evolving in a spatially flat Friedmann-Robertson- Walker 
(FEW) universe with a background fluid which has an 
equation of state ps — wbPb- There exists just two 
possible late time attractor solutions with quite different 
properties, depending on the values of A and wb '■ 

(1) A 2 > 3(wb + 1). The late time attractor is 
one where the scalar field mimics the evolution of the 
barotropic fluid with wq — wb, and the relation Qq = 
3(w B + 1)/A 2 holds. 

(2) A 2 < 3(u>b + 1). The late time attractor is the scalar 
field dominated solution (£Iq — 1) with wq = — l + A 2 /3. 

Given that single exponential terms can lead to one of 
the above scaling solutions, then it should follow that a 
combination of the above regimes should allow for a sce- 
nario where the universe can evolve through a radiation- 
matter regime (attractor 1) and at some recent epoch 
evolve into the scalar field dominated regime (attractor 
2). We will show that this does in fact occur for a wide 
range of initial conditions. To provide a concrete example 
consider the following potential for a scalar field Q: 



V(Q) =M 4 (e aKy + e pKW ) 



^^Q\ 



(1) 



where we assume a to be positive (the case a < can 
always be obtained taking Q — > — Q). We also require 
a > 5.5, a constraint coming from the nucleosynthesis 
bounds on Qq mentioned earlier |l7|Jl§| ]. 

First, we assume that (3 is also positive. In order to 
have an idea of what the value of (3 should be, note that 
if today we were in the regime dominated by the scalar 
field (i.e. attractor 2), then in order to satisfy observa- 
tional constraints for the quintessence equation of state 
(i.e. wq < —0.8), we must have (3 < 0.8. We are not 
obviously in the dominant regime today but in the tran- 
sition between the two regimes so this is just a central 
value to be considered. In Fig. [I] we show that accept- 
able solutions to Einstein's equations in the presence of 
radiation, matter and the quintessence field can be ac- 
commodated for a large range of parameters (a, (3). 

The value of M in Eq. (0) is chosen so that today pq ~ 
p c « 10~ 47 GeV 4 . This then implies M « 10- 31 A<f P1 « 
10 _3 eV. However, note that if we generalize the potential 
in Eq. (@) to 
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V(Q) = Ml, x {e aK(Q - A) +e pR(Q - B) ) , 



(2) 



FIG. 1. 
Contour plot of wq (today) as a function of (a, (3), with 
the constraint £Iq (today) ~ 0.7. The region a < 
5.5 is excluded because of the nucleosynthesis bound, 
f2(g(lMeV) < 0.13, and the upper region due to la ob- 
servational constraints. 



then all the parameters become of the order of the Planck 
scale. Since the scaling regime of exponential potentials 
does not depend upon its mass scale [i.e. M in Eq. (Q)], 
A is actually a free parameter that can, for simplicity, be 
set to Mpi or even to zero. On the other hand, just like 
before for M, B needs to be such that today we obtain 
the right value of pq. In other words, we require M 4 ~ 
M^e-f 33 - p Q . This turns out to be B = O(100)M P1 , 
depending on the precise values of a, (3 and A. 

There is another important advantage to the poten- 
tials of the form in Eq. (|l|) or Eq. (||) ; namely, we obtain 
acceptable solutions for a wider range of initial energy 
densities of the quintessence field than we would with say 
the inverse power law potentials. For example, in Fig. 
we show that it is perfectly acceptable to start with the 
energy density of the quintessence field above that of ra- 
diation, and still enter into a subdominant scaling regime 
at later times; however, this is an impossible feature in 
the context of inverse power law type potentials |2£| . 

Another manifestation of this wider class of solutions 
can be seen by considering the case where the field evo- 
lution began at the end of an initial period of inflation. 
In that case, as discussed in Ref. pq] , we could expect 
that the energy density of the system is equally divided 
among all the thousands of degrees of freedom in the cos- 
mological fluid. This equipartition of energy would imply 
that just after inflation fi^ w 10~ 3 . If this were the case, 
for inverse power law potentials, the power could not be 
smaller than 5 if the field was to reach the attractor by 
matter domination. Otherwise, Q would freeze at some 
value and simply act as a cosmological constant until the 
present (a perfectly acceptable scenario of course, but 
not as interesting). Such a bound on the power implies 
wq > —0.44 for fiQ — 0.7. With an exponential term, 
this constraint is considerably weakened. Using the fact 



that the field is frozen at a value Qf rs Ck — y/GTTl/K, 
where Qi is the initial value of the field [E5| , we can see 
that the equivalent problem only arises when 



2 In a > In 



PQi 

2Peq 



(3) 



where pQ i is the initial energy density of the scalar field 
and p eq is the background energy density at radiation- 
matter equality. For instance, for our plots with a,i = 



10- 



l eq 



= 10 4 , this results in a bound a < 10 3 . 




log a 
FIG. 2. 
Plot of the energy density, pq, for a — 20, (3 = 0.5 
and several initial conditions admitting an Qq = 0.7 flat 
universe today. The solid line represents the evolution 
which emerges from equipartition at the end of inflation 
and the dotted line represents p ma ttor + ^radiation- 

A new feature arises when we consider potentials of 
the form given in Eq. (|l|) with the nucleosynthesis bound 
a > 5.5 but taking this time (3 < 0. In this case the 
potential has a minimum at KQ m i n = hx(—/3/a)/(a — 0) 
with a corresponding value Kiin = Af 4 ^(-§) Q /( Q "/ 3 ). 

Far from the minimum, the scalar field scales as de- 
scribed above (attractor 1). However, when the field 
reaches the minimum, the effective cosmological constant 
Fmi n will quickly take over the evolution as the oscilla- 
tions are damped, driving the equation of state towards 
wq = — 1 . This scenario is illustrated in Fig. 0, where the 
evolution of the equation of state is shown and compared 
to the previous case with j3 > 0. In many ways this is the 
key result of the paper, as in this figure it is clearly seen 
that the field scales the radiation (w = 1/3) and mat- 
ter (w = 0) evolutions before settling in an accelerating 
(w < 0) expansion. Once again, as a result of the scal- 
ing behavior of attractor 1, it is clear that there exists 
a wide range of initial conditions that provide realistic 
results. The feature resembles the recent suggestions of 
Albrecht and Skordis pQ]. The same mechanism can be 
used to stabilize the dilaton in string theories where the 
minimum of the potential is fine-tuned to be zero rather 
than the non-zero value it has in these models Ib3| . 




log a 
FIG. 3. 
The late time evolution of the equation of state for pa- 
rameters (a, /3) : dashed line (20,0.5); solid line (20, —20) 
and flQ fa 0.7. (a = 1 today). 



In p5J, a quantity T = V"V/(V') 2 is proposed as an 
indicator of how well a given model converges to a tracker 
solution. If it remains nearly constant, then the solutions 
can converge to a tracker solution. It is easy to see from 
Eq. (n]) that apart from the transient regime where the 
solution evolves from attractor 1 to attractor 2, T = 1 to 
a high degree of accuracy. 

It is important to note that for this mechanism to work, 
we are not limited to potentials containing only two ex- 
ponential terms and one field. Indeed, all we require of 
the dynamics is to enter one period like regime 1, which 
can either be followed by one regime like 2, or by the field 
settling in a minimum with a non-zero vacuum energy. 
We can consider as an example the case of a potential 
depending on two fields of the form 



V(Qi,( 



M 



4f aiK,Qi+a 2 K,Q2 



/3ikQi+/3 2 k(; 



(4) 



where all the coeficients are positive. This leads to similar 
results to Eq. (§]) for a single field Q, with effective early 



and late slopes given by a^ s = 



and/3 2 ff = #+# 



'2) 

respectively Such a result is not surprising and is caused 
by the assisted behavior that can occur for multiple fields 
P5| . Note that for this type of multiple field examples 
the effective slopes in the resulting effective potential are 
larger than the individual slopes, a useful feature since 
we require a e s to be large. 



III. DISCUSSION 

So far, we have presented a series of potentials that 
can lead to the type of quintessence behavior capable 
of explaining the current data arising from high redshift 
type la supernovas, CMB and cluster measurements. The 
beautiful property of exponential potentials is that they 



lead to scaling solutions which can either mimic the back- 
ground fluid or dominate the background dynamics de- 
pending on the slope of the potential. We have used 
this to develop a picture where by simply combining po- 
tentials of different slopes, it is easy to obtain solutions 
which first enter a period of scaling through the radiation 
and matter domination eras and then smoothly evolve to 
dominate the energy density today. We have been able 
to demonstrate that the quintessence behavior occurred 
for a wide range of initial conditions of the field, whether 
p Q be initially higher or lower than p matt or + Eradiation- 
We have also shown that the favored observational val- 
ues for the equation of state wq (today) < —0.8 can 
be easily reached for natural values of the parameters 
in the potential. This is a big improvement in respect 
to most quintessence models as they usually give either 



[5] I. Zlatev, L. Wang and P.J. Steinhardt, Phys. Rev. Lett. 
82, 896 (1999). 



WQ 



> 



-0.8 or wq 



1. 



We have to ask, how sensible are such potentials? Can 
they be found in nature and, if so, can we make use of 
them here? The answer to the first question seems to 
be, yes they do arise in realistic particle physics models 
p6| [OJ , but the current models do not have the correct 
slopes. Unfortunately, the tight constraint emerging from 
nucleosynthesis, namely a > 5.5, is difficult to satisfy in 
the models considered to date which generally have a < 
1. It remains a challenge to see if such potentials with 
the required slopes can arise out of particle physics. One 
possibility is that the desirable slopes will be obtained 
from the assisted behavior when several fields are present 
as mentioned above. 

It is encouraging that the quintessence behavior re- 
quired to match current observations occurs for such sim- 
ple potentials. 
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